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In modern turbofan engines, variable geometry has been 
incorporated to improve some off -design performance. Most 
control designs ignore this variable geometry and use fuel 
metering as the primary control input. 

This thesis investigates the use of variable geometry to 
control the engine and, thereby, reduce fuel consumption due 
to transients. Additionally, steady-state trim conditions 
are altered to reduce the static fuel consumption. The non- 
linear transient simulation program is used to analyze the 
steady-state operating condition and develop small perturba- 
tion control limitations. Linear models, both large and 
reduced order, are used in analyzing the effect of controllers 
on system response. A computer program was generated to re- 
duce a large order linear model to a usable size for control 
system development. 

This analysis shows the reduced-order regime dependent 
controllers to be viable and to favorably enhance the quest 
for reducing specific fuel consumption in existing engines . 
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INTRODUCTION 



I . 

High performance, military aircraft are operationally 
tasked to conduct missions throughout their designed per- 
formance envelopes. To accomplish these missions, the air- 
craft must operate efficiently within its flight regime. 

The typical jet engine is designed for operation at a single 
specified altitude and Mach number combination, e.g., 30,000 
feet at M = 0.9. This design point will give inherently 
good cruise performance and is particularly well-suited to 
commercial transport operations. This engine design method 
forces the user into accepting off-design performance in 
typical military applications. Engine performance needs to 
be optimized to efficiently conduct all assigned missions, 
e.g., air combat maneuvering, ground attack, supersonic inter- 
ception and carrier operations. Adaptive digital control 
methods can be incorporated within the engine to accomplish 
this mission dependent optimization. 

Modem control methods allow the designer to develop a 
control system to regulate every parameter of the parent 
system. He has been aided by computer simulations of non- 
linear systems and linear approximations of those systems. 

The desired control system can be developed parallel to or 
ahead of the controlled system. 

The state-of-the-art turbofan engine utilizes variable 
geometry, in such forms as variable incidence inlet guide 
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vanes, adjustable exhaust nozzle area and airflow bleed-off. 
Hydromechanical controls have been used in the past to im- 
prove some elements of engine performance through positioning 
these variable engine components. Electronic controls, such 
as were introduced in the Pratt & Whitney F100-PW-100 turbo- 
fan engine, allow improved scheduling of the variable devices, 
again improving the engine performance. Studies have been 
done on the FI 00 engine to optimize its performance at a 
point other than its design point [1] . 

Originally, the onboard microprocessing capability of an 
aircraft was intended for weapon system improvement and 
integration. Further computing equipment was added as part 
of the engine control to supervise the hydromechanical con- 
trol system and detect faults. The F100 engine introduced an 
electronic control which allowed use of digital control tech- 
niques incorporating the hydromechanical devices as implemen- 
tors rather than controllers. Digital device improvements 
continue to expand the available computing power available, 
making flight control and engine control integration a dis- 
tinct possibility. As more software space becomes available 
to the engine controller, regime adaptive control will be the 
accepted standard in military applications. 

This thesis looks at using modern, multivariate control 
logic in designing flight-regime dependent control for a 
typical military turbofan engine, the Pratt & Whitney F100. 
Specifically, the attempt is made to improve thrust specific 
fuel consumption, i.e., fuel economy, in the high altitude 
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cruise environment. Such improvement in fuel economy is 
especially important in the conservation-conscious decades 
ahead. 
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II. TURBOFAN ENGINE PERFORMANCE CONSIDERATIONS 



A. FLIGHT REGIME FACTORS 

The ability of a military aircraft to perform an assigned 
mission is determined by its flight regime. Figure II-l is 
a typical flight regime for a modern, high-performance 
military aircraft. 

The confines of the flight regime are dictated by numer- 
ous aerodynamic and propulsive factors. A well designed 
aircraft will incorporate an engine having operating boundar- 
ies in excess of those of the airframe. This design technique 
allows for a factor of safety in actual operation as well as 
compensate for installation and equipment aging losses. 

At the low speed, low altitude region of the envelope, 
engine performance is judged in the ability to produce 
sufficient thrust for take-off and landing and to provide 
rapid engine acceleration in the missed approach or balked 
landing. The extreme upper region of the flight envelope 
requires the engine to produce maximum thrust while not ex- 
ceeding turbine temperature limitations. 

Figure II-2 is a mission delineated flight regime for the 
supersonic tactical fighter aircraft. The desired mission 
performance from this chart must be interrelated with the 
aerodynamic and propulsion limitations of Figure II-l. An 
adaptive control system can be designed to schedule optimum 
performance of the airframe and engine throughout the flight 
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Figure Il-l. Flight Envelope Limitations 
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Figure II-2. Flight Regimes 
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regime. Northrup Corporation is involved with a U.S. Air 
Force project to implement an adaptive flight control system 
in an aircraft such as the F-16. The same type of approach 
can be applied to engine control and integrated into an 
aircraft control system to optimize the total package. 

B. TURBOFAN ENGINE PERFORMANCE MEASUREMENTS 

The primary factor with which to judge engine performance 
is the engine's ability to convert fuel to thrust efficiently. 
This measurement, thrust specific fuel consumption, SFC, is 
a function of the airflow through the engine, M , the bypass 

cl 

ratio, b, and the fan and core exhaust velocities, U ' ' and 

e 

U ', also termed cold and hot stream velocities, 
e 



where M^ is the fuel flow, f is the fuel-to-air ratio and V 
is the flight velocity. Figure II-3 shows the F100-PW-100 
engine in cross section with reference labels. 

The engine designer fixes the airflow and bypass ratio of 
the engine. The exhaust velocities are functions of the ex- 
haust temperatures and pressures, which are controllable 
through fuel metering and variable geometry. 



M 



SFC 



M^ (1+f)U 
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Figure II-3. Engine Cutaway 
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where y and y ^ are the respective ratios of specific heats, 

T q6 and T Q8 the respective exhaust total temperatures, n is 
the exhaust nozzle efficiency, and P a /Pgg and P a /Pg8 
respective ratios of free stream ambient pressure to dis- 
charge total pressure. 

The standard bypass fan will produce a fixed pressure 
ratio and exhaust temperature, producing an exhaust velocity 
in the cold stream in excess of the flight velocity. This 
increases the engine's thrust without altering the hot stream 
through the engine core (ignoring energy lost in driving the 
fan) . The advanced, low-bypass fans that reintroduce the 
cold flow into the augmentor, as in the F100, increase the 
augmentor efficiency and during non-augmented operation, 
allow for increasing U e ' 1 through the use of the variable ex- 
haust nozzle. 

Studies show that Thrust Specific Fuel Consumption, SFC, 

can be decreased through increased turbine inlet temperature 

and compressor pressure ratio. Both of these variables are 

designed into the engine at the engine design point. Engine 

control hardware can be designed to maintain the turbine inlet 

temperature, FTIT, and pressure ratio, R , at the optimum or 

c 

best achievable level during off-design operations. 

C. ENGINE CONTROL DEVICES 

The mechanical devices available for engine control are 
the hydromechanical fuel control, the variable exhaust nozzle 
and variable incidence inlet guide vanes or stators. The 
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first device has the effect of increasing the fuel-to-air 
ratio and temperatures in the turbine and exhaust sections of 
the engine. The combined effect of increasing f and U ' is 
increased* thrust , but does not necessarily improve SFC due to 
the corresponding increase in the fuel flow. The variable 
area exhaust nozzle has been used for the last two decades in 
most high performance military aircraft. This device in- 
creases the hot section exhaust velocity, U 0 ' , thus increasing 
thrust and maintaining near ideal pressure in the exhaust 
section of the engine. 

The latest mechanical devices are the variable incidence 
blading concepts. The first use of these devices was to 
achieve near design mass flow rate through the compressor 
during the engine start cycle. These variable stators would 
be rotated to a starting configuration until the engine reached 
a given rotational speed and then be driven to a fixed normal 
operating position. 

The full advantage of variable geometry blading comes in 
the optimizing of off -design performance. The single stage 
pressure ratio is found to be 

p n , U 2 C 

= 1 + n pf— [ 1 - A tan b. + tan a.)] (II-4) 

P 01 RT 01 r 2 1 

where n is the stage efficiency, U is the rotational velocity 
of the blading, R is the Universal Gas Constant and Tq^ the 
total temperature at the front face of the stage. At the 
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design condition, the factor (tan b 2 + tan a^) is minimized, 
thus providing the best pressure ratio. The angle b 2 is the 
relative flow direction upon leaving the rotor portion of the 
machine; a^ is the angle at which the flow impinges the rotor 
and is a function of the stator angle (see Fig. II-4) . Con- 
trolling the stator angle increases the stage pressure ratio 
in the early stages and increases the compressor (or fan) 
performance over a wide range of operations. 

In the F100-PW-100 engine, four mechanical devices are 
used, fuel control, variable exhaust nozzle, inlet guide vanes 
with moveable trailing edge and moveable stators in the first 
three compressor stages. The exhaust nozzle maintains nozzle 
area, expansion ratio and boattail drag, simultaneously, near 
optimum. The inlet guide vanes are used ahead of the fan to 
improve inlet distortion tolerance, improve fan efficiency 
and enhance engine acceleration performance. The variable 
stators improve starting and high Mach number characteristics . 

D. ENGINE MODELING 

It is not economically practical to build an engine and 
then conduct experimentation to determine its performance. 
Modern technology, in both digital and hybrid computers, 
allows the design engineer to build a simulation model of the 
engine and predict the system response well in advance of 
hardware manufacture. Several generic computer simulations 
are in use, these allow non-specific engine performance analy- 
sis. Additionally, a computer model is generated for each 
new development engine by the engine manufacturer. 
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II-4. Compressor Stage Flow 
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Szuch, [2] , provides a detailed description of current 
engine modeling techniques. The Pratt & Whitney F100-PW-100 
engine has been modeled with both a real-time hybrid computer 
[3] , and a digital transient simulation [4] . Both programs 
have distinct advantages and applications. Although the digi- 
tal program is not real-time, it provides a large volume of 
data and is easily accessed and maintained. A copy of the 
F100-PW-100 Transient Engine Simulation Deck (CCD 1103.3.0) 
was obtained from the YF100 Special Project Office of the 
Aeronautical Systems Division of the U.S. Air Force at Wright- 
Patterson Air Force Base, Ohio. This program was modified 
to be compatible with the Naval Postgraduate School IBM 
370/3033 computer system. 

Designing a controller based on the non-linear digital 
simulation alone is extremely difficult, if not impossible. 
Miller and Hackney [5] , developed high order linear approxi- 
mations, at steady state operating points, from the non- 
linear simulation program. The linear system was modeled in 
the classical state variable form 

X = FX + GU (II-5) 

Y = HX + DU ( II-6) 

The 16 system states, listed in Table II-l, were chosen to 
coincide with the states measured in the non-linear simulation. 
These states include 2 rotational speeds, 3 internal total 
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TABLE II-l 



xl - 
x2 - 
x3 - 
x4 - 
x5 - 
x6 

x7 - 
x8 

x9 - 
xlO - 
xll - 
xl2 - 

xl3 - 

xl4 - 
xl5 - 
x!6 — 



Engine States 

Fan speed, (Nl) , in RPM 
Compressor speed, (N2) , in RPM 
Compressor discharge pressure, (Pt3) , in psia 
Interturbine volume pressure, (Pt4.5), in psia 
Augmentor pressure, (Pt7m) , in psia 

Fan inside diameter discharge temperature, (Tt2.5h), 
in Rankine 

Duct temperature, (Tt2.5c), in R 

Compressor discharge temperature, (Tt3) , in R 

Burner exit fast response temperature, (Tt4hi) , in R 

Burner exit slow response temperature, (Tt41o) , in R 

Burnerexit total temperature, (Tt4) , in R 

Fan turbine inlet fast response temperature, (Tt4.5hi), 
in R 

Fan turbine inlet slow response temperature, (Tt4.51o), 
in R 

Fan turbine exit temperature, (Tt5) , in R 
Duct exit temperature, (Tt6c) , in R 
Duct exit temperature, (Tt7m) , in R 
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pressures and 11 total temperatures. The inputs to the sys- 
tem are the fuel flow, the 3 variable geometry elements pre- 
viously discussed, and the customer bleed air percentage. 

The outputs include the net thrust, the engine total airflow, 
turbine inlet temperature and the stall margins. Addition- 
ally, two fan exit pressure ratios are determined for some 
of the operating points. The inputs and outputs are listed 
in Table II-2. 

The technique used by Miller and Hackney perturbed each 
state, xi, slightly while holding the other states and inputs 
constant. This allows calculation of the deviations caused 
by that state in the other states and outputs. These devia- 
tions then form the F and H matrices of equations (II-5) 
and (II-6) . Perturbing each input, ui , in the steady state 
configuration determine the G and D matrices. The sampling 
time was set at 7 milliseconds. 

The linear models obtained are very good approximations 
of the non-linear system but are often too complex and/or do 
not contain the most convenient parameterization to be used 
in a practical design. A reduction of model order and, per- 
haps, augmentation of the linear model is generally required 
prior to using the model for control design. 

E. CONTROL DEVICE EFFECTIVENESS 

Before designing an engine control, one must determine 
the effect each input has on the system and the responsive- 
ness of each input. The five control inputs listed in Table 
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TABLE II-2 



ul - 

u2 

u3 

u4 

u5 

yl - 
y2 - 

y3 - 
y4 - 
y5 - 
y6 - 

y? - 



Engine Inputs 

Main burner fuel flow, lb/hr 

Nozzle jet area, ft**2 

Inlet guide vane position, deg 

High compressor variable vane position, deg 

Customer bleed flow, % 

Engine Outputs 
Net thrust, lb 

Total engine airflow, Ib/sec 

Turbine inlet temperature, R 

Fan stall margin 

Compressor stall margin 

Fan exit delta p ratio (test data) 

Fan exit delta p ratio (theory) 
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II-2 are used in controlling the F100 engine. The customer 
bleed air percentage, u5, is essentially constant in steady 
state operations and can therefore be ignored in the control 
development. A study could be made to determine the merit 
in extracting additional bleed air from the engine as a pri- 
mary control, and dumping that portion not required by the 
customer systems. 

The intent of this control development is to minimize the 
specific fuel consumption of the engine. At the same time, 
only small variations in turbine inlet temperature, y3, and 
total engine airflow, y2, are desirable and no change in 
thrust is allowable. One method might employ using exhaust 
nozzle area, u2 , and fan inlet guide vanes, u3, inputs to 
minimize SFC at the trim condition and then compensate for 
thrust loss by using fuel to balance the equation. 

Figure II-5 plots the change in specific fuel consumption 
due to incremental changes shown in Table II- 3 in the inputs 
ul, u2 , and u3; input u4 , the compressor variable guide vanes, 
exhibit no influence on SFC. Proportional derivatives, i.e., 
slopes were determined by making a linear approximation to 
the curves. It is clear in the figure that nozzle area and 
inlet guide vane inputs (u2 and u3) have the most substantial 
effect on SFC. 

Figure II-6 plots the net engine thrust, yl, versus the 
control input changes. Again u2 has the greatest effect and 
ul and u3 contribute about equally to thrust variations. It 
is also noted that a linear combination of ul and u3 can be 
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n-5 . Change in SiC vs Contxrol Inputs 
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Figure II-6. Change in Thrust vs Control Inputs 
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Control Inputs 



TABLE II -3 
Input Perturbations 

The following perturbations are used throughout this 
investigation. They represent the largest expected value, 
or physical limitation placed on each control input. 



Fuel Flow, ul 
Nozzle Area, u2 
Inlet Guide Vanes, u3 
Compressor Vanes, u4 
Bleed Air, u5 



± 200 lb/hr 
± 0.2 sq . f t . 
± 4 degrees 
± 4 degrees 
0.0 percent 



used to offset the thrust change caused by a change in u2 . 

This is important in meeting the zero net thrust change de- 
sign requirement. Figure II-7 shows the effect the inputs 
have on turbine inlet temperature, y3. Again, u2 is the 
dominant control input. These figures clearly demonstrate 
the need for minimizing the exhaust pressure mismatch with 
ambient conditions. 

Table II-4 is a summary of the proportional derivatives 
extracted from the preceding figures. Using these derivatives 
the following relationships are found. 

d(SFC) = 2.27E-05 d(ul) + .0655d(u2) 

- . 00092d (u3) ( II-7 ) 

d (y 1) = .9325 d(ul) - 1351.7d(u2) 

+ 18 . 6275d (u3) (II-8) 
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TABLE II -4 

Summary of Derivatives 



Input 


SFC 


Thrust (yl) 


FTIT (y 3) 


Fuel Flow, ul 


2.27E-5 


.9325 


.12126 


Nozzle Area, u2 


.0655 


-1351.7 


283.085 


Inlet Guide Vanes, u3 


-.00092 


18.6275 


-2.1641 



U4 was found to have no effect on the parameters of 
interest . 



d(y3) = . 12126d (ul) + 283.085d(u2) - 2.1641d(u3) (II-9) 

If a decrease in u2 of -.2 square feet and an increase of +4 
degrees in u3 are applied, a change of -.017 occurs in SFC. 

The fuel flow, ul, is used to balance the thrust equation. 

To accomplish this the fuel is reduced by 370 pounds per hour, 
further decreasing specific fuel consumption. The total effect 
of these three inputs at the trim condition is a 3.5% reduc- 
tion in SFC, no thrust change and a 110 degree decrease in 
the turbine inlet temperature. 

By requiring that thrust level not change, the SFC per- 
formance is compared under identical operating conditions. 

Each input can be studied separately to determine its effect 
on the outputs and on SFC. To remain in the small perturba- 
tion regime inputs should not be varied more than 10% of their 
steady-state values. Combinations of inputs are then chosen 
to effect the most desirable improvement in SFC and minimize 
the effect on engine performance. 
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This technique has shown improvement in the specific fuel 
consumption at the operating point chosen in the studies of 
Reference 1. The use of the existing variable geometry to 
improve the fuel efficiency at the high altitude cruise 
operating point chosen for this investigation will be dis- 
cussed in a later section. 
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III. CONTROL DEVELOPMENT 



A. LINEAR OPTIMAL CONTROL THEORY 



Given the linear system, 



X 



FX + GU 



(III-l) 



consider a linear control law defined by 



U 



-KX + U 



0 



(HI-2) 



The original system can be modified to 



X 



(F - GK) X + GU q 



(III-3) 



The response of the modified system, as well as its eigen- 
values, is determined by the effect of the gain matrix, K. 

For a low order system, the control gains could be found by 
trial-and-error methods and applied to the system to determine 
if the response was acceptable. This method is impractical 
for large order systems. 

Optimal control is achieved by minimizing a cost function, 
J, which is defined as 



where R^ and R 2 are symmetric weighing matrices; R^ is non- 
negative definite and R 2 is positive definite. A third portion 
can be added to ( I II —4 ) to provide terminal state weighing. 



J 



J / [ (X T R 1 X) + (U T R 2 U) ] dt 



( III-4) 
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(III-5) 



x x (t 1 )P 1 x(t 1 ) 



where is a non-negative definite symmetric matrix. 

It is shown by Kwakernaak and Sivan [8], that the optimal 
input, Ug(t), based on the cost function is given by 



U Q (t) = (R 2 ) -1 G T P(t) 



(III-6) 



thus the knowledge of p(t) solves the regulator problem. A 
2n linear system is formed with the optimal system behavior, 
Xg(t), and the adjoint variable, p(t). 



x Q (t) 



p(t) 



-G r 2 G- 



• R 1 ' F ' 



n 


x 0 (t) 




p(t) 

— J 



(III-7) 



where 



p(t) = P (t ) X Q ( t) 



(III-3) 



with 



K(t) = R 2 1 (t)G T (t)P(t) (III-9) 

P(t) is a symmetric non-negative definite matrix that 
satisfies the Riccati equation 
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-p(t) 



R x (t) - P ( t ) G ( t ) R~ 1 ( t ) G^ ( t ) P ( t ) 
+ P ( t) F (t ) + F T (t)P(t) 



(III-10) 



with the terminal condition of 




P 



1 



(III-ll) 



Under these conditions, as the terminal time, tl, approaches 
infinity, the control law asymptotically approaches a unique, 
stable steady state condition. Thus, the matrix Riccati 
equation, (III-10) , can be rewritten 



Many numerical methods are available to solve the matrix 
Riccati equation for the steady-state control gains. Some 
of these are covered by Kwakernaak and Sivan [6] . This 
thesis employed a computer solution, 0PTSYS4 , developed by 
Hall and Bryson [7], at Stanford University. A representative 
output is included in Appendix 3. 

B. CONTROL SYSTEM MODELS 

At this point, an analysis of control models is made to 
establish the form of the controller to be applied. A 



0 



-IT T 

R -PGR. G P + PF + F P 

1 2 



(III-12) 



and 



K 




(III-13) 
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comparison of small and large perturbation controls must be 
made. Take the linear system 

X = FX + GU (III-14) 

Y = HX + DU (III-15) 

This system can be made closed-loop by incorporation of a 
control law 



U = - KX . ( III-16) 

K can be a time-varying matrix or constant as required by 
the system. A more general control law would be 



U = U + C (X - X ) + / C A(Y - Yjdt (HI-17) 

m x m J y u 

where A is the diagonal matrix whose diagonal terms are 1 
or 0 determining which outputs are used in formulating the 
input. and C are again matrices associated with the 

states and the outputs, respectively. Additionally define 



U 



m 



U n + R t 
0 u 



(III-18) 



X 



m 



X. + R t 
0 x 



(III-19) 



the 0 subscript indicates trim position/value; and R x are 
functions of u and x associated with time. For small pertur- 
bations the time varying terms, R u and R , can be ignored and 
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equation (III-17) can be rewritten 

6 U = C 6X + / C A6Y dt (III-20) 

x ; y 

Dropping the delta terminology and taking the derivative of 
this and using equation (III-16) , yields 

• • 

- KX = C X + C AY ( III -21) 

x y 

combining terms 

(-K-C ) X = C AY ( III-22 ) 

x y 

Substituting, using the system description of equations 
(III-14) and (III-15) one obtains 

(-K-C ) (FX+GU) = C A (HX+DU) (III-23) 

x y 

and solving for C^, 

C = -K-C A(H-DK) (F-GK) -1 (III-24) 

x y 

and using all outputs (A = I) and equally weighing them 
(Cy = I) one has 

C = -K- (H-DK) (F-GK) _1 (III-25) 

X 
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consider- 



However, if the outputs are ignored (A = 0), i.e., 
ing the case where small changes are made in the output 
vector , then 



as in the simplified control law of equation (III-16) . 

This approach does not apply to the large perturbation 
case. In large perturbation analysis, the time varying 
terms cannot be ignored. The mean values of u and x are as 
in equations (III-18) and (III-19) . and are the time 

functions that describe the path that the mean values follow 
over the period of the large perturbation. The control law 
becomes 



Again, if all outputs are used and equally weighed, then 
taking the derivative 



C 



x 



K 



(III— 26) 



5U -R t 
u 



C (SX -R t) + 1C AoY dt 
x x ; y 



(III-27) 




CX-CR -CRt+CAY 
x xx xx y 



(III-28) 



and combining terms 



-C (F-GK)X + R C +R C t 

X XX XX 



K (F-GK) X + R, 



u 



+ R t + (H-DK) X 



(III-29) 
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If we set 




R 



u 



(III-30) 



and 




R. 



u 



(III-31) 



in effect, assuming a linear relationship between the states 
and inputs . Then 



as in equation (III-25) for the small perturbation case. 

An adaptive control system removes the need to model the 
large perturbation response. Linear models can be developed, 
as described in Section II. E, to define the large perturba- 
tion case as a series of linear relations. The onboard digi- 
tal computer then schedules the control gains based on the 
linear models as the flight conditions pass through the model 
transition points. The rapidity of response inherent to the 
digital computer removes the requirement for modeling the 
total dynamic system, rather, only the dominant response model 
for the selected operating point. The reduced-order linear 
model is far easier to handle in formulating the control sys- 
tem and is derived to match the dominant system characteristics. 



-C 



- K - (H-DK) (F-GK) 



-1 



(III-32) 



x 
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C. REDUCTION OF LINEAR MODELS 



DeHoff and Hall [8], present a method of reducing the 
order of the linear system based on the dominance of the 
states and the crosscoupling of the states as observed through 
eigenvector analysis. The required order of reduction is 
determined by the number of dominant states appearing in the 
desired control bandwidth. For the F100 study, DeHoff and 
Hall chose a control bandwidth of 0 to 10 rad/sec, corres- 
ponding to the primary control device (fuel flow actuator) . 
Their operating point (static, sea-level, intermediate power) 
dictated that the control be designed to modulate thrust, 
providing maximum thrust for take-off performance. At another 
operating point a different parameter may dictate the control 
design and, thus, the primary control actuator and control 
bandwidth. 

The method of decomposition involves finding the eigen- 
vector transform matrix, T, such that 

XT = TA (III-33) 

where A is the diagonal eigenvalue matrix. Then, by defining 
an alternate state vector, Z, where 

X = TZ (III-34) 

the original system may be written in modal coordinates as 

Z = AZ + T~ 1 EU (III-35) 
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If a control bandwidth is determined, the eigenvalues 
falling in that range are the only ones that need to be in- 
cluded in the initial reduced order model. The matrix may 
be reordered by constructing a matrix Q, which has a 1 in 
the position corresponding to the eigenvalue to be included 
(column) and the reordered position (row) of that eigenvalue. 
Thus, if eigenvalue 5 of the original A matrix was to be re- 
ordered to the first position, a 1 would be placed in the 
(1,5) position of the Q matrix. Now, the reordered matrix 

A' = QA (III-36) 

is incorporated in the linear model 

QZ = A'Z + QT _1 BU (III-37) 



letting 



Z ' = QZ 



(III-38) 



Z ' 



A ' Q _1 Z ’ 



QT~ 1 BU 



(III-39) 



and 



X = TQ -1 Z' 

A matrix R can be found such that 



(III-40) 



X' = RX = TZ' 



X 



1 



X 



2 



J 



( III— 41) 
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This X' vector is the reordered state vector, where the 
elements contained in the subvector X.^ are those states 
associated with the eigenvalues in the control bandwidth 
and X 2 contains all other states . 

The linear system may be rewritten 



X' = RFR _1 X ' + RGU 
Y = HR _1 X ' + DU 

Now define 

F ' = RFR -1 

G ' = RG 

H' = HR -1 

and 




I 

H’ = [H. i H_] 

1 l 2 . 



(III-42) 

(III-43) 

(III-44) 

(III-45) 

(III-46) 



(III-47) 



(III-48) 



(III-49) 
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If the states in truly model the full order system, 
then X 2 will be essentially zero and 



X 1 = F 11 X 1 + F 12 X 2 + G 1 U 



(III-50) 



0 ’ F 21 X 1 + F 22 X 2 + G 2 U 



(III-51) 



Solving for X in equation III-51, one has 



X 2 - - F 22 1f 21 X 1 + F 22 1g 2 U 



(III-52) 



This is incorporated into the first equation 



X 1 - (F 11- F 12 F 22 1f 21 )X 1 + (G l- F 12 F 2i lG 2» U 



(III-53) 



and the output equation becomes 



Y = (H 1- H 2 P 22 1f 21 )X 1 + <D - H 2 F 22 1 G 2 , ° 



(III-54) 



The reduced order model can be written as 



X = FrX x + GrU 



(III-55) 



Y = HrX n + DrU 



(III-56 ) 



where 



Fr F H " F 12 F 22 F 21 



(III-57) 
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Gr 




(III-58) 



Hr 




(III-59) 



Dr 




(III-60) 



This technique was used to develop the program REDUCR2 , 
documented in Appendix A. REDUCR2 was used to develop the 
reduced order models used in the control development portion 
of this thesis. Sample input and output of REDUCR2 are 
listed in Appendix A. 

D. REDUCED ORDER MODELS 

Three reduced order models were generated with the REDUCR2 
computer program (Appendix A) . The criteria for selecting 
the desired states was that used by DeHoff and Hall, a con- 
trol bandwidth of 0-10 Hertz. The eigenvalues of the full 
16th order system were found using the OPTSYS4 computer pro- 
gram and then associated to the system states using modal 
analysis. Table III-l gives the results of this modal analy- 
sis. The five states associated with the control bandwidth 
are fan speed (xl) , compressor speed (x2) , augmentor pressure 
(x5) , burner exit temperature (xlO) , and fan turbine inlet 
temperature (xl3) . 

A fifth order model was generated using these states. It 
was noted that only three of these states (the two rotational 
speeds and the pressure) were dynamically significant, i.e.. 
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TABLE III-l 



Operating Point Eigenvalues 



State 


Physical Property 


Eigenvalue 


X4 


P04.5 


-540.23 


X3 


P03 


-156.432 


Xll 


T04 


-49.69 ± j 13 . 99 


X16 


T07 




X9 


T04Hi 


- 54.12 


X12 


T04 . 5Hi 


- 42.61 


X7 


T02 . 5C 


- 22.59 


X14 


T05 


-21.58 ± j 2 . 51 


X15 


T06C 




X6 


T02.5H 


- 18.62 


X8 


T03 


- 14.49 


X5 


P07 


- 7.35 


XI 


Nl 


3.60 control 

band- 


XI 3 


T04 . 5LO 


1.98 width 


X2 


N2 


- 1.31 


X10 


N04LO 


. 647 
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TABLE III -2 



Reduced Model Comparisons 





Model #1 


Model #3 


Model #2 


Full Ord 


XI 


173.51 


175.89 


175.89 


173.41 


X2 


125.10 


126.86 


126.86 


125.02 


X5 


.326 


.325 


.325 


. 326 


X10 


4.85 






4.86 


Xll 




53.51 




53.48 


XI 3 


3.98 






3.98 


XI 6 




23.04 




23.03 


Yl 


94.48 


174.8 


97.70 


156.99 


Y3 


11.84 


53.51 


9.3 


53.48 
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represent energy storage in the form of torque and pressure. 

A third order model was generated using these three states 
and a comparison made. In this initial analysis these two 
models produced virtually identical response with respect 
to the key outputs thrust (yl) and turbine inlet temperature 
(y3), see Table III-2. The 0PTSYS4 program was used to deter- 
mine the transfer functions associated with these two outputs 
for each model. It was found that the fifth order model con- 
tained zeros near the poles at -1.960 and -.651, the eigen- 
values associated with the two temperatures in the fifth 
order model. This caused these two poles to have no effect 
on the response of the model. Thus the third order system 
sufficiently models the control bandwidth. 

In both models, the turbine inlet temperature (y3) response 
did not closely model that of the full order system. There- 
fore, a third model was formed using the three states of 
the third order model plus the conjugate pair of poles repre- 
senting the turbine inlet temperature (xll) and exit tempera- 
ture (xl6) . This model gave excellent thrust and temperature 
response (yl and y3 , Table III-2) due mainly to the inclusion 
of the turbine inlet temperature as a system state. Table 
III-3 lists the three models, their states and associated 
eigenvalues . 

Although the conjugate pair of eigenvalues lies outside 
the control bandwidth, their inclusion markedly improves the 
output correspondence with respect to the full order system, 
see Table III-2. Further analysis will use the third order 
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TABLE III-3 



Reduced Model Eigenvalues 





Model #1 


Model #3 


Model #2 


States 


XI 


XI 


XI 




X 2 


X2 


X2 




X5 


X5 


X5 




X10 


Xll 






X13 


XI 6 




Eigen- 


-7.175 


-51.65 ± 


-6.947 


values 


-3.392 


j 7 . 14 


-3.282 




-1.418 


- 1.388 


-1.389 




- .651 


- 6.909 






-1.960 


- 3.265 
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TABLE I II -4 



Fifth Order Model 
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TABLE III -5 



Third Order Model 
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model (model #2) and the fifth order augmented model (model 
#3) . These models are listed in Tables III-4 and III-5. 

E. VERIFICATION OF REDUCED ORDER MODELS 

The reduced order model must be sufficiently representa- 
tive of the original system for the designed control to pro- 
vide the desired full system response. Weinberg and Adams 
[9] , used the F100 non-linear simulation to develop a linear 
17th order system at the static, sea-level point at partial 
power settings. Their reductions to 5th and 3rd order and 
similar reduction to 5th order by DeHoff [10], from Miller 
and Hackney's 16th order model show excellent correlation to 
the non-linear simulation. 

The operating point chosen for consideration in this 
thesis is 30,000 feet, Mach = 0.9 and power level angle of 
67 degrees. This corresponds to a near optimum cruise point 
and is a condition frequently encountered in routine opera- 
tions. Section III.D lists the states and associated eigen- 
values at this operating point and discussed several reduced 
order models. The dominant eigenvalues were considered to 
be those which affected the models fast response and had 
dynamic significance to the engine's physical characteristics. 
Again, this implies a control bandwidth of 0-10 rad/sec, how- 
ever, other models, using eigenvalues outside the control 
bandwidth, were considered in determining the best reduced 
model for deriving the control laws. 

One method to analyze model correspondence is to use the 
Bode plot across the control bandwidth. Figure III-l is a 
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Figure III-l. Model Comparison in Frequency Domain 



Bode comparison of the DeHoff reduced model and the Miller 
and Hackney 16th order model at static, sea-level, inter- 
mediate power for the 0 to 10 rad/sec control bandwidth. 

It is clear that through the 0-10 Hertz bandwidth the two 
models show very close correspondence; this similarity breaks 
down in the higher frequency region, but the response time of 
the eigenvalues outside the control bandwidth is considered 
to be too short to affect transient behavior. The phase 
response, not shown, follows the same trends. 
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IV. ANALYSIS OF CONTROL SYSTEMS 



A. COST FUNCTION WEIGHING FACTORS 



The control law for the engine regulating system, 



U 




(IV-1) 



can be reduced as shown in Section III.B to 



U 



U. 



m 



KX 



(IV- 2) 



when considering small perturbations and ignoring the out- 
puts from the system. In most circumstances, the U m , the 
trim condition, would be considered as zero for steady state 
analysis. The linear models developed by Miller and Hackney 
[5], set all inputs to zero at the steady-state operating point 
for which the model is defined. 

The K matrix is the solution to the optimum regulator 
problem discussed in Section III. A. The 0PTSYS4 computer 
solution used in this thesis uses the cost function 



replaces the need for the R^ matrix in equation (III-7) . 

Kirk [11] , discusses the effect that weighting matrices 
have on system performance. In general, the larger the 



J 



\ f [Y T A y Y + U T R 2 U]dt 



(IV-3) 
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magnitude of the matrix element the faster its associated 
parameter (system state, input or output) stabilizes. An 
additional factor in determining the magnitude of the diagonal 
matrix elements is the magnitude of their associated param- 
eters. Given the cost function 



J 



1 

2 



/ X T R 1 X dt 



(IV-4) 



where 



II 


1 

a 2 


0 


; x = 


xl' 




0 


1 

b 2 




x2 






D 







The cost function then becomes 

J = j l(~) 2 + dt ( IV- 6) 



The desire is to minimize the cost function, J, with 
respect to time. Thus, taking the time derivative and 
setting it equal to zero 

x 2 y9 2 

(~) + ra = o (IV-7 ) 

a b 

This is the solution which produces a minimum cost. The 
values of a and b are the largest expected values of xl and 
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x2, respectively, and represent an equal effect on the cost 
function by each state. Should either a or b be decreased, 
the optimal control developed would have to increase the 
control input to the other state to minimize the cost func- 
tion. A nominal fluctuation limit of ±5% will be considered 
as a small perturbation in the initial control analysis. 

Designing a control to minimize the thrust specific fuel 
consumption (SFC) of the engine requires limiting the thrust 
and the fuel flow fluctuations. The weighting matrix ele- 
ments corresponding to these two parameters are, therefore, 
to be increased with respect to the other parameters. Looking 
at only the small perturbation response, it is noted that 
very little change in the total engine airflow ( y 3 ) is ex- 
pected and, as before, there is no change in customer bleed 
air requirements (u5) . These two parameters can be neglected, 
as their response characteristics have no effect on the desired 
SFC response. To demonstrate the effect of the different ele- 
ment magnitudes on the control gain matrix, the 0PTSYS4 pro- 
gram was run for the 3 cases delineated in Table IV-1. Case 
I was developed to allow each output and input, with the excep- 
tion of the two pressure ratios (y6 and y7) and the bleed air 
percentage (u5) , to effect equal contribution to the cost 
function. The diagonal elements of both A and R ? normalize 
the squared value to the order of one. Case II uses all of 
the same elements, except the thrust (yl) and the fuel flow 
(ul) are weighted more heavily, to the magnitude of ten. This 
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TABLE IV- 1 
Case Definition 



A and R 
Y 

elements are 



are diagonal matrices, only the significant 
listed . 



Case I 



A 

y 


= 


( .0001 


1 .0001 100 100 le-8 le-8) 


R 2 


= 


( . 0001 


100 .01 .01 le-8) 








Case II 


A 

y 


= 


(.01 1 


.0001 100 100 le-8 le-8) 


R 2 


rr 


(.01 100 .01 .01 le-8 ) 








Case III 


A y 


= 


(1 1 


.0001 100 100 le-8 le-8) 


R 2 


= 


(1 100 


.01 .01 le-8) 



heavier weighting should reduce the fuel input to the system 
and moderate the thrust fluctuations as the control returns 
the system to the steady state condition. Case III increases 
the weighting of yl and ul to achieve a normalized magnitude 
of 100. This extremely heavy weighting should eliminate fuel 
from the control and reduce thrust variations to absolute 
minimum. The whole intention of the controller design is to 
damp out state perturbations without using fuel or altering 
the thrust level, i.e., making no change in specific fuel 
consumption to effect control. 
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These three cases were implemented on the two reduced 
order models and simulations run for a perturbation of -50 
RPM in fan speed (xl) . The fifth order model returned the 
fan speed to 95% of steady state in .5 seconds using Case I, 
but increased fuel flow by 10 pounds per hour initially and 
required an increase in fuel consumption of .2% over the 
half second period. For Case II, the fifth order model re- 
turned to 95% of the steady state fan speed in approximately 
.3 seconds, while using almost 80% less fuel than the Case I 
controller. This 80% reduction in fuel consumption indicates 
that the increased weighting of Case II has accomplished 
exactly what was desired. Case III was run for the fifth 
order model and, although virtually no fuel was used in the 
control, the system experienced oscillations in fan speed of 
±8 RPM through the 3 second simulation run. Additionally, 
the Case III controller caused rapid changes of ±4 degrees in 
the compressor guide vanes (u4) , the regulatory limit imposed 
for the simulation, and also large fluctuations in the inlet 
guide vanes (u3) . One would achieve results of similar magni- 
tude for each of the cases assuming a +1% perturbation of all 
states in the reduced-order model. Case I is too fuel de- 
pendent and causes large changes in u3 , Case III also requires 
large changes in U3 that could cause a slew rate request larger 
than that available from the control actuator. 

In modeling the control, the actuator and physical limi- 
tations of each control must be considered and appropriate 
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limitations placed on the control law to avoid overshoots. 
Table IV-2 gives the actuator limits and maximum rates 
applicable to the F100 engine [4]. A stepping limit of .01 
seconds is assumed in determining the rate limits and fuel 
flow fluctuations are limited to ±10% of the steady-state 
value. No compensation is allowed for hysteresis or non- 
linear dynamics of the actuators; these are considered 
minimal in the small perturbation case, but would need to 
be incorporated for large perturbation control modeling. 



Input 
Fuel flow 
Nozzle area 
Inlet guide vanes 
Compressor vanes 



TABLE IV-2 

Actuator Limitations 

Maximum Minimum 

16300 lb/hr 450 lb/hr 

6.4 sq. ft. 2.8 sq. ft. 



0 deg. -40 deg. 

4 deg. -40 deg. 



Rate Limit 
15800 lb/hr/sec 
3.6 sq.ft. /sec 
48 deg/sec 
40 deg/sec 



It is evident that the Case III weighting is not accepta- 
ble because of the oscillations and the rapid, large amplitude 
fluctuations in the guide vane positions (u3 and u4) . Case I 
is also inappropriate for the specific fuel consumption mini- 
mizing control since it uses fuel as the primary control in- 
put. For these reasons, Case II is chosen as the weighting 
for all control modeling to be implemented in later sections. 

Once the integral portion of equation (IV-1) and the feed- 
back gains are determined, the constant portion, U , must be 
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analyzed. In most situations this term would be neglected. 

The steady state analysis is now performed to determine if 
the specific fuel consumption can be reduced by making small 
changes to the trim position of the variable geometry. An 
analysis of control changes versus SFC, as was done at the 
static sea-level condition in Section II. D, must again be 
conducted. A minimum SFC obtainable relates to certain input 
values in U m which is incorporated into the control law, im- 
proving the steady-state performance and controlling the 
fluctuations of the system which will have the greatest effect 
on the SFC. 

As was done at the static sea-level operating point, a 
small perturbation in each input was made and its effect on 
the SFC, thrust and other engine conditions tabulated. The 
goal of these analyses is to develop a trim input, U , that 
improves the fuel economy at steady state without causing 
significant change in the operating conditions (both states 
and outputs) . 

The effects of the control inputs on several parameters 
are plotted on Figures IV-1 through IV-4. Figures IV-1 
shows, again, the large effect nozzle area (u2) and fuel flow 
(ul) have on specific fuel consumption, the effect of inlet 
guide vanes (u3) is drastically reduced from that observed 
at static, sea level operation. Figure IV-2 indicates the 
offsetting effect of ul and u2 and again indicates a zero 
net thrust control can be found. Fan speed (xl) and compressor 
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Figure IV-1. SFC vs Control Inputs 
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Figure IV-2. Thrust vs Control Inputs 
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Figure IV-3. Fan Speed vs Control Inputs 
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Figure IV-4. Compressor Speeds vs Control Inputs 



speed (x2) effects are shown in Figures IV- 3 and IV- 4, respec- 
tively. At this operating point the inlet guide vane posi- 
tion (u3) may be most useful in controlling the fan speed 
fluctuation and thereby stabilizing the fan stability margin 
(y4) . A sample control would close u2 by .2 square feet and 
reduce fuel flow by 57 pounds per hour. This combination 
produces a 1.4% decrease in specific fuel consumption and 
results in a 4% reduction in fan speed. A trim condition 
control input, U m , would be 
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B. REDUCED ORDER SIMULATION ANALYSIS 

The model developed in Section III.C and Section III.D for 
the high altitude flight operating point was used to initially 
analyze the control law developed in Section IV. A. A reduced 
order controller is simpler and cheaper to implement through 
a digital computer. It requires less computer space and far 
less computing time than the full-order linear model developed 
controller. 

A comparison of the three cases for which gain matrices 
were found indicated that Case 2 gave the most desirable re- 
sults. The Case 2 weighting matrix was applied to the third 
order reduced system and perturbations applied to various 
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Figure IV-5. Third Order versus Fifth Order 
Control Response 



states. The same weighting matrix was applied to the fifth 
order model (#3) and similar perturbations applied. Figure 
IV-5 shows the system response to a 50 RPM perturbation of 
the fan speed. The fifth order system has a settling time 
of approximately .45 seconds whereas the third order system 
has one of 1.5 seconds. This variation is accounted for by 
the two fast response poles in the fifth order model, 

(-51. 65, j7. 14) . 

These results must be compared to the 16th order system 
to determine which model most closely approximates the full 
order system response. 

C. FULL ORDER SIMULATION ANALYSIS 

The proof of the viability of a reduced order controller 
is its ability to control the full order system in the same 
manner as predicted in the reduced order model. The two 
reduced models analyzed in Section IV. B were applied to the 
sixteenth order linear model. Additionally, the Case 2 weight- 
ing matrix was applied and a full order controller was found 
and implemented. All simulation was done using the IBM Con- 
tinuous System Modelling Program, Model III. 

All three models were compared with similar perturbations 
applied. Figure IV-6 compares the full order and fifth order 
response to the 50 RPM speed perturbation. Figure IV-7 com- 
pares the full order and third order response to the same 
perturbation. Figure IV-8 compares the two reduced order 
controllers and the full order controller for the same fan 
speed perturbation. 
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Figure IV-6. Fifth Order vs Full Order Control Response 
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Figure IV-7. Third Order vs Full Order Control Response 
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Figure IV-8. Full Order Implementation 



The results show satisfactory agreement between the re- 
duced order and full order implementation of the reduced 
order controller. This proves the viability of the con- 
troller and implies that comparable results will be found 
if the control is implemented in the non-linear dynamic 
system. The fifth order controller, with two poles outside 
the control bandwidth, most closely models the response of 
the full order controller. As was seen previously, the 
inclusion of the conjugate pair of poles dramatically im- 
proves the response time of the model and provides a much 
closer approximation to the full order system response. 

This indicates that the fifth order model is sufficient 
for the control development process at this operating point. 
The two fast response poles outside the control bandwidth 
must be included to achieve this close simulation. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The study of implementing a reduced order controller 
that is flight regime dependent shows that at a specific 
operating point an improvement in specific fuel consumption 
can be made. The conclusions reached at the end of this 
research are: 

1. The variable geometry existing in the state-of-the- 
art turbofan engines can be used to further improve the 
specific fuel consumption; 

2. A reduced order controller can be found and imple- 
mented that concentrates control authority on the variable 
geometry rather than the fuel control, thus furthering the 
fuel savings; and 

3. A similar controller can be found and implemented at 
any given operating point throughout the flight regime, these 
controllers could then be organized into a single digital 
machine to provide a regime-dependent engine control to im- 
prove the airframe/engine interaction in all flight conditions. 

B. RECOMMENDATIONS FOR FUTURE STUDY 

One obvious extension of this work is to find further 
reduced order controllers. 

Additionally, work can be done to fully implement these 
reduced order controllers in the non-linear model using 
either the Pratt & Whitney transient simulation digital model 
or the hybrid computer model. 
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The effect on fatigue life and life cycle cost due to 
the "fuel saving" control schemes presented was not deter- 
mined. It is therefore recommended that these factors be 
explored in order to fully analyze the usefulness of these 
control schemes. 
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READ IN THE FULL ORDER MODEL MATRICES AND THE EIGENVECTOR MATRIX 
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MODAL MEASUREMENT SCALING MATRIX 
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THE CONTROL WEIGHTING MATRIX 



CO 














O 














1 

Q 














O 














o 














OOOOo 














• • ♦ • • 














OOOO—* 














CM 














o 














1 

Q 












^ N ro xf --t 


o 












ooooo 


o 


• 


• 








(fill 


o ooo^ 


• 


• 








cocon 


• • • • • 


X 


oo 








— <r-vr oco 


000^0 


UJ 


UJ •• 








inmtnrocg 




►— 


oo 




CM— ILT\ 




cneneg— irn 




m 


.JO 




ooo 




— 4 com I s - 




>• 


<T + 




1 1 1 




O rr O— 1 — ' 


CM 


00 


>o 




QQO 






o 




20 




aoco— i 




• • • • • 


1 


CL 


UJ—- < 




e- vQC7> 




—4 — H —-4 r—4 


O 


o 


O — » 




enj-'O 




1 1 


o 


a 


—•m 




moo va- 






o 




ujcr 




in cr %o 




^CO'OOO 


OOOOO 




• 




• • • 




O—IOOO 


• • • • • 


Q 


oC cm 




in err- 




1111 + 


oo^oo 


OJ 


O 1 




1 


• 


ooooo 




<S) 


>- *• 






• 


4*-4roieO 




o 


<tm 


• 


— icn 


X 


m co<f^o 




_J 


JO 


• 


ooo 


•— « 


CMO00— *00 


CM 


o 


o + 


X 


i i i 


ac 


>00-100 


o 




OQ 


t—l 


ooo 


V- 


(MO'stO'O 


+ 


-J 


UJ 


QC 


a — i>o 




njOO— *— • 


o 


< 


ocr^ 


L~ 


rg<rc0 




• • • • • 


o 


X 


CM 


< 


mmO 




mmro^o 


o 


»~4 


-JO 


X 


CMOLT, 


a: 


1 1 1 


ooooo 


h- 


<1—4 




vOQN 


a 






CL 


X • 


ct 


• • • 


h- 


o — <r-e-cM 


OHOOC 


o 




o 


co in eg 


o 


o— lOOO 






| 


h- 


1 1 


UJ 


\ 1 1 1 + 




UL 


a «• 


o 




> 


ooooo 




o 


OvO 


MJ 


— 4— < eg 


X 


oomcno 


M- 




o 


> 


ooo 


UJ 


iT\N sj* >0 eg 


o 


X 


CL + 




1 f f 


O 


cOvtOOn 


1 


UJ 


oo 


IU 


OOO 


►— « 


sTLOO^Jf^- 


a 


h- 


oo 


O 


LOsOO 


UJ 


«-* sO vj-tn cm 


o 


L/? 


-JO 


— * 


'to sO 




nOGOCM^T o 


o 


> 


CO 


UJ 


inmr^ 


-J 








OCM 




njoom 


C 


O vOCM^J— * 


• • • • • 


z: 


UJO 


►— 


vr 


QC 


1 1 1 1 


-•o^oo 


jj 


<S) • 


X 


• • • 


1— 






o 


o<n 


o 


oo '“■■* 


X 






— « 


-J 1 




1 f ! 


c 






UJ 


o 


nc 




o 





97 



CLDSEO LOOP OPT. REG. LEFT EIGENVECTOR MATRIX 
- I . 840600-0 l 5. 364860-03 -5. 66l 900+01 
- I . 012120*0.) 6 . 0489 20-02 5. 738020+01 

-o. 032820-01 L. 035080+00 1.408710+00 
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CSMPIII SAMPLE PROGRAM 



in 

•» 

r- 

Q 



o 






— 














•k 






sO ^ ^ ^ 














— k 




O 


— •r— NOtn 














in 




• 


« fr 














•» 




o 


*— 4 ~4 •> *>— 4 














nO 




«k 


|| || ^4*-4 t| 














•—4 




o 


— II II »“• 














w 




o • 


c 

I 

} 

1 














CQ 




O i O'? 
















•k 


— * 


o • •» 


vO vO — — vO 






o 








— -k 


in 


• oo 


h -4 in in —» 


«■» 




cn 








sO 


w 


O sO • 
















0—4 


V— 


* g-rg 


r4 k— 4 *—4 ■— 4»— 4 


k— 4 


—4 


a 








m 


•• 


co oo in g- - 


M ll il li ii 


— 


O' 


i— 








o 


— k 


1 l — • **rg 


—>—> -)~i 


> 












—4 


— rg 


OC *Oo<^ • 




4 


O 


o 








— 1 


f— — 


OO »QQ • * 


— « — k — «k 


^«k 


— 1- 


o 


— * 


— k 




< 


—a: 


• • O * *(\| • 




rg 






“5 


X 


“5 


•» 


X •* 


hhi^Oo^O 


1 

1 

< 

< 

< 


— 


—a 


— k 


— 


— 




a 


*^-k 


* * Oin 


M 


at: 


xo 


— k 


X 


o 


X 


vO 


— r- 


ii it ii tnr-g-^ 




— 


# 


— k 


* 


4 — * 


■4 






<oocoa 






*— • 


—* 


*k.i. i 


a 


- • — > 


— il ii ii ii 


—*W 


k 




"k— 


o 




— > 



ntno •* 'OvOO 



— O— ■ >in*:u. i i | margin 

O— 1)11 

♦ O— ••aixxw^-^ 

ct — u.uj — xo >a:>- 

UJLO *0 Oo^ct 
o 2 — <QOct:ct UJUJIUUJLJU 
UJ UJ O OC UJ x UJ tJU _J _J _J _J 
>- x — »c x h- .jm coco co coco 

UJLUCLX<<I<< 
J— O LLimU. srac <1— f— H-»— H- 



in— - 



^ooooo 

UJOOOOO 
O eg rg rg rgrg 

i— in ininmin 

IQQQQC 
<X< <<<< 
QUJUJUJUJUJ 
U-a: Q^QCQiCt o 



—4 | •— « >M ^ 

— -Owl— — — >0 >C — — — 

oin •— «w «—4 —«<inco— 

4 •» ..4- UJ— UJH- •. ..4 «k4 — 

-4 ii o ii — • ii • ii uj ir co ii — mi — • ii 

— • O — — ■ » — —»».-» 3 — • • O— — — ■ 

CC O CO UJ IL~ 

— o It m II — — O II II o II O II ^ 



o — 

— <OlTv 



O — • — 

U.O— C'-'U. LL 

-t-^OOOO'H — 



oc — -a— c 

UQUJU-CUJO 



— u 
—a 

ILV 



-•co 

UJ II OC 

x oo 

•“•O II II 

4 ^0 — k— k 

— — 

oc— 

ucoc 



II — II 

o 

O il o 

r- — ro 

0 — 0 
no^ 





cc 


3 " a: 














H- 


oo 


<o 














k— i 


mo 
















n** 


Org 


>a 


n 


-» c 


o 


o 


o 


«n> 


/ 

/ 

/ 

i 


2 


QZ 


cr 


CT' O 


g* 


n 


m 





99 



k. 



OOOO OOCOOOoOr OCOOCOCCO^COOC 
+ + + + 4- + + + +4-I | |4--L4.+ |4.4. + J4.| | \ ± 

OQQGQOOOCQQOOQOOCQQOOQCOOCO 
CCOOOO^OOCOOCOCCOOOOOO n OOrc 
uOO'Cgu\sOO>J-^nr< >m rx JO sOsC'O 4 — *0 uO'^C 
O^o^ovin^jv^rcsjcoac t o '.f cr v r ajcco'-o j- 
o^r— <\jxj- 0 s o Lrrginm.— o roOM'n--i o o ^ . 
O-f m— * — «r i^Occ >f jin— * — ‘^JT r J •r. 



OOOOOOOoOOOOOOOoO^OO^r: ooooo 

I I I I I I I I I 1 fill I i 

OOOOOOOOOCOOOOOOOCOOCOOCOCC 

+ + ♦+ + + + + + + I | |4.|4.+ |*.4-4.|+|+ + 4- 

QnQnoooQoooQOQoaGGoncionnnno 
OOOOOOOOOOOO^OO^OCCOOOOCOCO 
coo comr-^rsia > ocsjod^gq-nT u^oo > oro.j-m.oin^ro7 v ccr 

j''t , x\^ r co r o^Nf 

— » rHCMT.^lT. r fSJC*' cnjvlt.O cccroo Vocco 

O v— •r\i—» in r\i r\jc\j o — — <lt — «aj_4irun —^lt. — * 

o 

• OOOO ooococcocooc ooooc occo^c 

il I I I I I I II I I I I I I I 

_J OOOOOOOOOOOOOOOO oooooooocoo 

UJ + +I+ + + + + + + +I I | 4. I ♦ ^4. ^ * I + I I + 1 

a oooo oooo oooooooo oooo oocoono 

• OoOOOO OOOOOOO OOOOOOO ooooooo 

in oa'Crfxioo^ir»r-r-f\jvt tou^vroN^M-^cc^^ njir. m 

o K -o^oo s ^onj'Or-^OLnvO^^'f cr — ^ 

• Oin "OO X><N—*<r r^^i^o^ou^cr occc Qmv.>oco 

il <r<? o— » -Hrno^^^or-r- o— ccir^tMri^vij^OsO^ 



o 

u. 

on 



X 

3 



QX 

+ * 

«OJ 

— O 
102 
H II — 

— ~C 



a 

o 

X 



a: 

o 

H- 

X 

II 



-*2 




2T 


t: 


r-a 


o 


a 


a 


1 •— 1 




»— • 


•— t 


— < ►— 


1— 


►- 


V— 




< 


< 


<f 




-J 


-j 


-J 


-o 


D 


3 


3 




3 : 


s: 


3: 


in^irg 


M-4 


*— i 


►— « 


1 on—* 
— * X 


on 


on 


LO 


— LU ► 


UJ 


UJ 


LU 


02- 


o— . 


on 


U1 




mm* 




NH 


— O'— 


O'- 


3ll 


3 


o^x 


QC> 


a :on 


a: 


—4 UJ » 


o - 


O - 


o 


1 LU 1— UJH- UJ1— 



-<02003:00200 
^O m JO*^jO M JO 

X — «l— Q »V*CL — ♦►-a.— ‘ 

LL VU- >LL >LL 
H- h~X j-x »-X 
3UJ3 LUO UJ3 lu 
CL -J CL UJ -J a UJ -j CL UJ _J 



0^0 <*— 3<»—3<~ 
Cl— OCL>— CCL 1—CCLl— 



Tv*,) ^ 


L- 


t: 




a: 


taMri 


1— 


3 


X 


QC 


LO 


oc: 


o OC 


3 


LU 


ao Oon 







LU OOOOOOOOOOOOOOOO 000000000^0 

O I I I I I I I I I I I I 

cn — 1^-iOrO -40— <—*000—4— *OO^m 

CL OOOO OOOO OOOO OOOO OOOO OOOOOOO 

•- 4-* + ++ ++ + + + + +|4-+*+|J4.+ | | J | 4- | 

iTt OOOOOOnOOQOOOOOOCOCOOOCOCOO 

O OOOOOOOOOOOOOOOOOOOOOOOOOOO 

• aomin— * oia no<t t <"vio*- no om ooocnr*- om 

ii r*-o^GNa s c\ioao^vivi“0<jo oinrjmcosoa ocxior* — 'M-O't 

«j inrnrncDOro^(Mrsjtnf\j^^-i NOO' 4 jtr > -ocvjsOrr> 0 ^<J'Cr • 

lu ^in^rsjrgooLn— rsio * 4 “^ — ‘ltvo— 'OO 

o 

►— OOOOOOOOOOOOOOOOOOOOOOOOOOO 

3 111 I I I I I I I I I ! I 

a oooo^of\»o<N— <o— *rg— *oo— *— iooo^o— <o 

- OOOOOOOOOOOOOOOOOOOOOOOOOOO 

O ++4-+4.4»*+|+4-*t | J | + + | + + 4- J | | | | 

• OOOO OOOO OOOO OOOO OOOOOOO OOOO 

m OOOOOOOOOOOOOOOOOOOOOOOOOOO 

il oom^oo -o— rr ^ ^incoiNno- moo— »ro 

T ao— ‘p-sT 'Drntn^'004 n o — ‘-or- r o<~\i o o o <l- • >-o 

— 1 00 CM/> vt vt' or^OO r Mvf^OOOO'4-f\J oo lAvT^vf^^XN 

»— rgv^rOrsj-rr\jrnin o— •— *oo^ — ‘ojin— *rg<r 

z • • • 

— OOOOOOOOOOOOOOOOOOOOOOOOOOO 

UL III II III III I III II I 

— * O *-* rO *— < O Oi ^ o O «-* ^ r\j r j r j — * o O ^ J — 1 O ^ J <* J o <“• j n, J 
oc ^OOOOOOOOOOOOOOOOOOOOOOOOOOO 

LU 3+ + 4.4.4. + -L4.4-+4.I \ |^|^4.4.^|4-4-^| | | 

SQCLGoOOQOOQOOQOQOOCQOnCOOCCCOO 

_j< — • z: ^ o o o o o o c o o o o o o o o o o o o o ^ o o o o o o 

<h-LLi^-,Mr4(M4 , mir\f r lvOOrorgOO s r\j sO — ’ O o*-4" >0 LT — * •> f 

L^cO^OLrNO~rOvOro^— 4 , 00 — ‘t’— 
rsjsOcofM^T'-*— <NOOO( > J 4 , a N rn(T'- *— nr^j < ^'n — *c j 
r j m <j* — ♦ >o 0- ^ J o j ^ r\j — < — * vj- r-- tT' q ^ —t j 



o oooo ooooooo OOOO oo^ 
M I III l 



I I I 



s ooo 

I I I 



J 



100 



pg~-4m— 4— * 0 — *rgo~-«-^rgO —TMin vf sT 

oooooooo OOOOOOO ooooo 

4. + + +II++ + 4.+ II4.+ -*-1111 

qqooqqqqqqqqqqq qocqo 

OCOOOOOOOOOOOOO OOOOO 

— 4mr*-r^ — 4 00 0 <Mr- cgrgvO — *0 cc'aOO'i" 

•cooo-'o-tmcrooo— TO ONrOvOm 

m rsj <t" o — 4 >0 rg r\j <t" —4 rvj m nj co m r- i r — < —4 O * 



m<f-m<r vT cv m <r <r 
OOOOO OOOOO 
I I I I I I I I I I 
QQOQQQQOOQ 

oooooooo oo 

— <<r— i^uMncMvornsO 
r^acrc?Ovr'to^( r i 

*^ 0 ^'OO^OQO S (\J«H 

s O -J" mrg-f vf— 4 <Min— *— • 



m 


rgrgo 


Org— 4C gCn 


O 


OOo 


OOOOOOO 


1 


1 1 4- 


f ++++++ 


o 


OOO 


QQOOOQO 


o 


OOO 


OOOOOOO 


-4 


vrmo 


coo s <\if^ mr^m 


M 


m<ro 


o g-^inm-4 


—4 


on m 


OiTivOriin- *rg 




omrg vf 


m o rg — 4 — 4 m m 


OOOOO 


ooooooo 



OoOOO OOOOO ooooo 

i i i i 1 * I 

<M rg m —4 ^ -4 — »0 *-* <M<m m — i 

oocoo OOOOO OOOOO 

+ 4. + +|+++| + | I ♦ + + 
QQOOOOOQQOQQQQC 
OCOOOOOOOOOOOOO 
'O rg m o m O < 4 * oo — 4 ro 0 s 0 s r- ~4 
4*r*'Tvj^ s 0'-0— *— ♦^'maorjo-^ 1 ^ 
O *-< — * oo O <M 'O ^ fM -+ r- CT >0 oo 
-4 eg c* nO — * -* co ro ^ r- —• 



O OOOOO OOOOO OOOOOOO 
I II II I I I I I 

rg^rgg-m *-<Nt^mnrnvfm!Min 

OOOOO OOOOOOOOOO 

♦ + I I I I I I I I 1 I I I I 

OOOOO OOOOOOOOOO 

OOOOO OOOOOOOOOO 

cf o oor-igo -*<aomr*-o s ocr N 'Ococo 

-Oin vO v^r\| go co^ooq^o^nC^ O 

t"*’ O 4* — * ~4 mo°°^rg — 4 -4rgm vO 

^^vOrr^OOcnjcgrj^mmrO^C 



I 

O 

+ 

o 

o 

o 

o 

ro 



I I I 



I I 



(\jfNJr^On^nO-^^0 
oooo OOOOOOO 

I I I + + + + ♦ I ♦ + 

OOOOOOO QOOO 
OOOO OOOOOOO 
NQmmaOs*- ocsj— i O' 

NrnoO'f^ — * *+ m 4* 

*-* o <r rg n- m —• g- O' m 



cro o^gD— *cg-^— i— 400— 4^4 



OOOOO OOOOOOOOOO ooooo OOOO OOOOO oooo oo OOOO OOOOOOO 



I I I I I I till I I I It 
og^o^rgrgorO— im^rgmocgg-g- 
OOOOOOOO ooooo oooo oooo 

+ + + +I++I++I + I+ + + + +I I I 

000000000000000000000 
ooooo ooooo OOOOOOO oooo 
-OsO^ or^g-orgmo' o mo m-^Osog^-^mm 
— -coccrgminocr mo'e-r-r- omoomrgrr v f-m 
>o m m co ^4 go m co m r - -4 r*- rg ao o rg ^4 go rg 

cr rsj >r in -4 g- g- -4 m <M -4 rg gj g* rg m cn —* o o o 



II i i i i i i i till t 

<M m g* <m g* g* g- g* ro g- g- -4 eg > f rg m g* -4 ro O rg rj — • g“ 



IIIIIIIIIII++IIII+++++I+ 

oooo no ooooooo oo on ooooooo 
o ooooo ooooo oooo oo oo oo ooo 

CGOgorgmCN-mcgincr — *g* mo mrginmegco c\ir- rg 
mor4mrgej<n4 _ r0f^-— *0— 4 Og">mg“f'^acmcr m eg g- 

0^( r '^0!Nm'0- < O'OOOr^-cof^'Om^r f\ir^ — »mco 

— « rg g* eg m rn vo g- g* r- -4 ro — « r*. in m eg so o m fMin — • 



OOOO OOOOOOO OOOO OOOOOOO OOOOOOOO OOOOOOO OOOO OOOOOOO 



I I I I I I I III 

4-<M— ‘OO — •ocgmrgegrg— <rg— ‘O— »rgm 

ooooo oooooooo oooo oooo 

++I++++I++I+I++++++I I 

oooo ooooo OOOO oooo OOOO 
ooooo OOOOO OOOOOOO oooo 
•t<NNsO^<MO«H^inN^Mm^^ina'Goo<7' 
commma'og}*-40r-rgc>0— «mg}N-g)r-N04' 

O' nT CM — < O cn h- mh- CO fOrg ^ 4 “ ^ o moo 

— • rg fM ir» vt ro <t — < o rg in rg ^ r*- go O < 



I I I I I I I I I I I I 

vt-^^—»vfii^vfrnrOLnvf-^rg 

O OOoOO OOOOOOO 
I I I I I I I I I I I + I 

oooo a oooo oooo 

O OOOOO OOOOOOO 
COf^^OO'Of^-^-MtrvGOOfM 
moo— ♦^rvjcornO'Osr^^ 
'g'mm0'4"0'4'<x>*4-m'4''0cr 
' O *t* <m o m<\jh- mm-^ oominvf 



II III 

4* rvjmrgr\jir\— « rgrg— < 
OOOOOOOOOO 
| I | + 4> + I + + <&. 

a ooooo oooo 

OOOOOOOOOO 
^moo^sO^ oon 
ooo^o -r >o omoo 
m^og^r ' — »og mu^-0 

OGO'Om^G'jpvj— 



OOOOO oooooooo oooooooo OOOO' 



I I I I I I I I I II 

f\jcg^— <cg-^o Ocg^cgo^ *-*0 <Mo 0 ~-*m<f ^ 

ooooo ooooooooooooooooo 

+ 4.4.g.^|4.|>^l + l+ g-4*^-f+l I I 
OOOOQOQOOOOOOOOOOOOOOO 
OOOOO o^o OOOOO OOOOO OOOO 
-^f^-rg<f ( McO'^0 N v^oo— ^vj-^oo^vf^oomm^m 
^ ^ ^ ^ o o o o m vo o so <m <r o o -o rg in oo 
oo >o -r in o 0 s m o rg vf >o o o 0 s 0 s C7^ m o go o vt* 



III I I I I 

vfNfvt*rgm^m^tinsr<rrgo 
oooo oooo OOOOO 
I I I I I I I t I I I + ♦ 

oooo ooooo oooo 

OOOOO OOOO OOOO 
Oh-^OoOrgrg-^vOO'mo-m 
O O rg ^ in oo orgin mm — 4-^ 
imrggo-4-rgsO-^ 



Or^OoOm^OOOO 
"l III 

O^nim— 4og<rorg— 4 
OOOOOOOOOO 

4^1 I ^ 4* +> 4. 4i 4- 4. 

OOOOOOOOOO 
OOOOOOOOOO 
— 4 — *omo— O'- 4 
h-c^mmmo— 'ccvco 
nOO— » rg— *M'f — <rg— < 
O ^1" — 4 - 4 rgco m — 4 nj cDvO 



co— 4m^m-4>^mNj“ — 4mrgvj-<Nj^ vfoo— ^oonimm^rgmrg^-trk 

ooooo oooo OOOO ooooo oooo oooo oooo 00000000000000^000 



I I I I I I I I I I I I I 

rsj— 4^pgr\jrgrg-^rg— 4 ^— 4 — <0— ^mm 

OOOO OOOO OOOOOOOOOO oooo 

4. 4- 4. 4* 4. I |^4>^|4>l 4 4» 4* | 4>| | | | 

ooooo oooo oooo oooo ooooo 
ooooooo coo oooooooo oooo 
r- o vf -4 vO r* o ^ ^ —4 r, m m m co — * m — « 

mrOvfmmcoo— -m-^fsjrvjoomogo^T^.foOM 
m r- rj go co o m o -^m n}- rg o ^ m so r- rjm m —4 
-4 ^-4 r g m H m rg m -vt* m g^ rg -4 m r- m rg co o m 



II I I I I 

^g-^i^rgg^^g'ing'rG^O 
OOOOO OOOOOOOO^ 

+ I I I I l I I I I i I I + 
OOOO OOOOOOOOOO 

OOOOOOOOOO OOOO 

O— 4 rg mo — 4 m nT— • oomoom 

ri/>m»4in^OC'N^'OrG-4N 
0>f'Ocomr fc -— 4vfG0rj— <og-rg 
O— < ^ —4 Mm rg»-4 n^Ng- — <m ^jo 



I I I I I I 

OO^— 4rv jrgOMOrg 
OOOOOOOOOO 

| 4» | I 4. 4 * + 4 4 

OOOOOOOOOO 
OOOOOOOOOO 
m *g- r* r j m — 4 m o>r 
Goommcoo'^—'— hA 
vO td — 4 m rg — * o cr ro 
m rg rvjrg — 4 r j vr ^ rg— • 



OOOOOOOOOO OOOOO'' 
III I I I I I I I 



'CO OOOO oooooooo OOOOOOOO OOOOOOOOOO 
I I I I I III II II I III 



101 



oooooo 

I 4 4 4 4 4 

oooaoo 

oooooo 

XOO'O'i’N't 
o or- oon 

P— ^ >0— <00 vf 



oooooo ooooooo 
oooooo oooooo o 
44444.44444444 
OQQO 30 Q OOOOOO 

oooooo ooooooo 
oooooo OOOOfMOO 

OOOOOOOOOOIACO 
OOOOOO ooooooo 
OOOOOO f>000 ^‘Or 



OOOOOO 

I I I II 
— *cn— iC\J 

oooooo 

144414 

oooooo 

oooooo 

nOoq OroOrn 
T'O'M- 

nj r\i o oo 

ITN f- <cOr\J 00 



000000 0000-1-00 

I 

sTOOOOO 00<M0r\i00 

oooooo ooooooo 

I + + f + l + f I + + f i 

OOOOOOO oooooo 

oooooo OOOOOOO 
OOOOOO OOOO^Oo 
— 'OOOCO OO 
(^COC^OOClTiO^OC 

oooooo ooooooo 



oooooo 

(till 

<Nmo~ 

oooooo 

I +1 I + + 

oooooo 

oooooo 

QN oo— «X) C\j 
r\j— kTKOO to 
O 00 O' vf O' — » 
Psi-^COmOJfM 



00000000^0^00 

ooooom ooooooo 
oooooo ooooooo 

4 4 4 + 4 I 4444444 

OQOOOOOOOOOCO 

oooooo ooooooo 
oooooo ooooooo 

C^^OO'OOOO^ oo^ 
ooooo^r ooooooo 
oooooo ooooooo 



oooooo 
I I I I 
.•f rg— 4^0 m 

oooooo 

4 4 4 14 4 

oooooo 

oooooo 

(Mf*~ iru+o ro 
r— T' — (7 s 

-0<\J — I'tCT' — * 



00000^0000 OOO 

ooo mo-+ ooo^Ovtc 
oooooo ooooooo 

4.4-4* I 4- | 4 4 4 4 4 4 4 

OOOQOOOQOO OOO 
OOOoOO OOOOOOO 
ooC'toin ooo 0 o<r>o 
ocONOtnooocono 
OCO^ONOOCCO^O 
OCOrOOflOoOOOO^O 



'"'OOOOo OCO^CMOOOOOinC 

I I I I I 

O^O-'OO^^OOoOO 00^00000 
OO OOOOOOO OOOOOOO OOOOO 
4 + I 44-4-1 444-444-4 I 444-44 

000000000000000000000 
00 000 OOOOO 000 000 000 00 

43v0'00— I r^r\l COOOOOO OOOOOOOO 
(MnJ* iT> vf ON N J - 'T' r 'OOCO OCCT'^OO^O 
xf^O'Oinst^COOCCOOOCOOCO 
1 r>a r\ir^-urn— *000000 00^40 0000 



OOOOOOO— *00 0000 OnT OOOOO 

I I I I I 

OCJIA — * — * O r\i o 00 oO o <\j OnJ- 00 OO o 

00 000 000 00 OOOOOO 00000 

4 | 4|444|4444^ | 4144444. 

cqcqoqqqqqqcqqqOgqcqq 

OOOOOOOO C'OOOOOOOOOO*— 
-rr-fM^OrnfM-^— • OO ooo^ O go OOOOO o 
MT n -<vO^^<T'OCO r 'O^^NO rv ^ r '' - ''l cO 

o-* err- -rr- 00 0000 o-r 00000 z u 

nTco ^-*infNj(NJoooooocooo-iooocoMa o 

• • • 

OOOOOOONOOOOOv0 n -4OOOCOZ^Z 

II II t I UJtOUJ 



102 






LIST OF REFERENCES 



1. Sain, M.K., and others. Alternatives for Linear Multi- 
variable Control , NEC Inc., Chicago, 111., 1978. 

2. Szuch, J.R. , "Techniques for Jet Engine Systems Modeling," 
Control and Dynamics , Volume 14, 1978. 

3. Szuch, J.R. and Seldner, K. , "A Real-Time Simulation of 
the F100-PW-100 Turbofan Engine Using the Hybrid Com- 
puter," NASA Technical Memo X-3261, NASA Lewis Research 
Center, 1975. 

4. Pratt & Whitney Aircraft, Florida Research and Develop- 
ment Center, "F100-PW-100(3) Transient Engine Simulation 
Deck. User's Manual for Deck CCD 1103-3.0," FR-9476, 
January 1978. 

5. Miller, R.J. and Hackney, R.D., "F100 Multivariable Con- 
trol System Engine Models/Design Criteria," Air Force 
Aero Propulsion Laboratory, TR-76-74, 1976. 

6. Kwakernaak, H. and Sivan, R. , Linear Optimal Control 
Systems , John Wiley & Sons, New York, 1972. 

7. Hall, W.E. and Bryson, A.E., "OPTSYS4 Computer Program," 
Stanford University, 1979. 

8. DeHoff, R.L. and Hall, W.E., "State Space Techniques 

and Modeling for Control," Control and Dynamics, Vol. 14, 
1978. 

9. Weinberg, M.S. and Adams, G.R., "Low Order Linearized 
Models of Turbine Engines," Aeronautical Systems Division 
TR-75-24 , July 1975. 

10. DeHoff, R.L. and Skira, C.A., "A Practical Approach to 
Linear Model Analysis for Multivariable Turbine Engine 
Control Design," Alternatives for Linear Multivariable 
Control , NEC Inc., Chicago, 111., 1978. 

11. Kirk, D.E., Optimal Control Theory , Prentice Hall, 
Englewood Cliffs, N.J., 1970. 

12. Hill, P.G. and Peterson, C.R., Mechanics and Thermo- 
dynamics of Propulsion, Addison-Wesley , Reading, Mass., 
1965 . 



103 



13. Speckhart, F.H. and Green, W.L. , A Guide to Using CSMP , 
Prentice Hall, Englewood Cliffs, N.J., 1976. 

14. Bryson, A.E. and Ho, Y.C., Applied Optimal Control , 
Hemisphere Publishing Co., Washington, D.C., 1975. 



104 



INITIAL DISTRIBUTION LIST 



No. Copies 

1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, Virginia 22314 

2. Library, Code 0142 2 

Naval Postgraduate School 

Monterey, California 93940 

3. Department Chairman, Code 67 2 

Department of Aeronautics 

Naval Postgraduate School 
Monterey, California 93940 

4. Professor D.J. Collins, Code 67Co 1 

Department of Aeronautics 

Naval Postgraduate School 
Monterey, California 93940 

5. LT Barry Lawrence Dougherty, USN 1 

Box 73 

East Meredith, New York 13757 

6. Mr. Charles A. Gardner 1 

Pacific Missile Test Center 

Point Mugu, CA 93042 

7. Dr. Gerhard Heiche 1 

AIR-310 

Naval Air Systems Command 
Washington, DC 20361 



105 



f 




uobb 
c. 1 



30 t. 



Oougherty 
eff ir BCt ° n fuel 

-ei :er Z?° l 
the cost tl0nS in 



i 



Thesis 1955:5 

D666 Dougherty 

c.l Effect on fuel 

efficiency of para- 
meter variations in 
the cost function for 
multivariable control 
of a turbofan engine. 



thesD666 

Effect on fuel efficiency of parameter v 




3 2768 002 00624 9 
DUDLEY KNOX LIBRARY 



